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Abstract 

We study the giant graviton solution as the A0IS5 x S 5 background is heated up to finite 
temperature. The analysis employs the thermal brane probe technique based on the blackfold 
approach. We focus mainly on the thermal giant graviton corresponding to a thermal D3- 
brane probe wrapped on an S 3 moving on the S 5 of the background at finite temperature. We 
find several interesting new effects, including that the thermal giant graviton has a minimal 
possible value for the angular momentum and correspondingly also a minimal possible radius 
of the S* 3 . We compute the free energy of the thermal giant graviton in the low temperature 
regime, which potentially could be compared to that of a thermal state on the gauge theory 
side. Moreover, we analyze the space of solutions and stability of the thermal giant graviton 
and find that, in parallel with the extremal case, there are two available solutions for a given 
temperature and angular momentum, one stable and one unstable. In order to write down the 
equations of motion, action and conserved charges for the thermal giant graviton we present 
a slight generalization of the blackfold formalism for charged black branes. Finally, we also 
briefly consider the thermal giant graviton moving in the AdSs part. 
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1 Introduction 



Recently, a new method to study thermal brane probes in string/M-theory has been explored. 
This method consists in using the blackfold approach [1, 2] 1 in the context of string/M-theory 
branes [6, 7, 8]. In the blackfold approach one can describe the dynamics of a black brane 
wrapped on a submanifold of the background space-time in the probe approximation where 
the black brane is much thinner than the length scale of the submanifold. This method has 
been applied to the thermalized version of the Blon system for the D3-brane [8, 9], the gravity 
dual of the rectangular Wilson loop as described by an F-string ending on the boundary of 
AdSs x S 5 [10], the M2-M5 version of the Blon system [11], and has been used to find a 
number of new black holes in String/M-theory [6, 7]. 2 

In this paper we apply the method to giant gravitons [14, 15, 16]. The archetypical case of 
a giant graviton is that of a D3-brane wrapped on a three-sphere and with the center of mass 
moving along the equator of the five-sphere in the AdSs x S 5 background. This corresponds 
to a blown up version of a point particle graviton and is described using the Dirac-Born- Infold 
(DBI) action for the extremal D3-brane. Our goal is to analyze what happens to the giant 
graviton as one heats up the AdSs x background to non-zero temperature, requiring the 
brane probe to thermalize with the background. We dub the resulting brane probe a thermal 
giant graviton. 

In the dual gauge theory description of AdSs x the giant graviton moving along the 
equator of 5 5 with angular momentum J is dual to a gauge theory multi-trace operator O gg 
with R-charge J and conformal dimension A = J. Heating up the AdSs x S 5 background 
and with it the giant graviton brane probe corresponds then on the gauge theory side to the 
thermal state that results from the ensemble of operators that are fluctuations around O gg . 
This is true up to the temperature Thp of the Hawking-Page transition where an AdS black 
hole is formed. Thus, having a description of a thermal giant graviton will provide important 
insight into the strong coupling behavior of the gauge theory side at finite temperature. 

We find in this paper the following features for the thermal giant graviton corresponding 
to a black D3-brane probe wrapped on a three-sphere with center of mass moving along the 
equator of 5 5 in the AdSs x •S' 5 background. Firstly, we find for a given temperature that 
the family of solutions has J > J m in > 0. Instead for the extremal giant graviton one can 
take the limit J — > thus connecting to the point particle graviton. The thermal giant 
graviton, however, is forced to be blown up to a finite-size three-sphere. This is analogous 
to what happens in the thermal Blon case where the throat of the brane - made of a D3- 
brane wrapped on an S 2 - has a minimal possible radius [8, 9]. Secondly, we find a maximal 
temperature T max which provides the scale of the temperature of the solution. However, the 
probe approximation gives T max 3> Thp, hence requiring T < Thp means that we have small 
temperature for the thermal giant graviton T/T max <C 1. This is analogous to the case of the 
thermal rectangular Wilson loop of [10]. Thirdly, for the free energy we find the following 

1 See also Ref. [3] which discussed the first application to neutral black rings in asymptotically flat space. 
Reviews include [4] and [5] gives a more general derivation of the blackfold effective theory. 

2 Large classes of new neutral black objects in asymptotically flat and (A)dS backgrounds were found in [12] 
and [13] respectively. 
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This is a prediction at strong coupling which can potentially be compared to what one can 
compute on the gauge theory side. Finally, we note that the phase structure of the family of 
solutions at finite temperature is similar to the extremal case in that for a given temperature 
and for J m \ n < J < J max one has two available solutions, one stable and one unstable. This 
is different than for the thermal Blon and the thermal rectangular Wilson loop cases in that 
there the number of available phases could become higher or lower when turning on the 
temperature [8, 9, 10]. 

Since we are employing the blackfold approach for studying the thermal giant graviton 
our black D3-brane probe consists of A^ coincident black D3-branes in the supergravity 
approximation with 7Vd3 3> 1 and N <C AA^. At the same time the probe approximation 
requires A^d3 "C N. 

In order to describe the thermal giant graviton using the blackfold approach we need to 
generalize the blackfold approach slightly. Firstly, we write down the blackfold equations of 
motion for charged branes in backgrounds with fluxes since the giant graviton corresponds 
to a black D3-brane in the AdSs x S 5 background in particular with a five-form RR-flux 
turned on. Secondly, we are describing a brane probe that moves with constant velocity along 
a Killing direction. Thus, it is not a stationary solution. However, the brane probe is not 
accelerating in that it is moving along a geodesic. Thus, we consider what we call quasi- 
stationary blackfolds in the sense that they correspond to boosted stationary blackfolds. In 
particular we describe the conserved energy and the momentum associated with the motion 
and we also extend the other physical quantities as well as the variational principle to the 
quasi-stationary case. We remark that since a quasi-stationary blackfold is not accelerating 
it does not emit radiation and one can thus go beyond the probe approximation and perform 
a matched asymptotic expansion for the full system of the background with the brane. 

For the extremal giant graviton there are two solutions for each value of J in the range 
< J < 9N/8 and the ones with the biggest S 13 radius are the stable solutions that minimize 
the energy (here we set the radii of AdSs and S 5 to one). Note that this family of solutions 
consists of two disconnected branches, one for < J < N and one for 1 < J < 9N/8. A 
similar feature is found for the thermal giant graviton in that for a given temperature there 
are two solutions for each value of J in the range J m \ n < J < J max and those with the biggest 
S 3 radius minimize the free energy and are thus thermodynamically stable. This is true both 
with respect to local and global stability. Thus we have shown that there are stable thermal 
giant gravitons in this range of J. Note also that the family of stable solutions again consists 
of two disconnected branches, one for J m \ a < J < N and one for iV < J < J max . 

In addition to thermalizing the "archetypical" giant graviton, i.e. a D3-brane wrapped on 
an moving inside the S 5 of AdSs x S 5 , one can also thermalize the other possible giant 
gravitons. In particular, we consider in this paper also briefly the D3-brane giant graviton 
moving on AdSs. The thermal effects observed in this case are similar to that of the S 5 case. 
There is again a temperature dependent lower bound on the minimum angular momentum, 
with no upper bound just as for the extremal case. Moreover, the small temperature expansion 
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of the stable branch of thermal giant gravitons in this case gives the same correction to the 
free energy as in (1.1). 

This paper is built up as follows. In Sec. 2 we review the extremal giant graviton configu- 
ration for a D3-brane expanded on the S 5 of the AdSs x S 5 background and set the notation 
for the rest of the paper. In addition to this, we also discuss a stable branch of solutions that 
has not received much attention in the literature. In Sec. 3 we discuss how to extend the 
blackfold method for branes in backgrounds with fluxes and derive the form of the conserved 
quantities and action for this case. We also introduce the notion of quasi-stationary black- 
folds. In Sec. 4 we find the thermal version of the giant graviton on S 5 using the blackfold 
approach and describe the solution space. We then subsequently make a detailed analysis of 
the thermodynamic properties of these thermal giant gravitons, including their stability, in 
Sec. 5. In particular, we derive the low temperature correction to the free energy given in 
(1.1). Finally, in Sec. 6 we consider very briefly the corresponding results for thermal giant 
gravitons on AdSs. We end with a conclusion and outlook in Sec. 7. 

A number of appendices is included providing further details. In App. A we present a 
detailed stability analysis of the extremal giant graviton solutions, including the non-BPS 
branch. The thermodynamical blackfold action for blackfolds in background fields with non- 
zero fluxes and a corresponding Smarr relation is discussed in App. B. App. C gives a derivation 
of the form of the blackfold equation of motion for thermal giant gravitons. Finally, in App. D 
we analyze the two meeting points of the thermal giant graviton solution branches. 

2 Giant graviton on S 5 revisited 

In this section we review the extremal giant graviton configuration in type IIB string theory 
on AdSs x S 5 . For definiteness we focus on the case in which the giant graviton sits on the 
S 5 , and correspondingly also the construction of the thermal version in Sec. 4 will be confined 
to this case. The case in which the graviton is expanded on the AdSs an d its thermal version 
is considered briefly in Sec. 6. 

The review below will serve to set our notation and properly define the configurations that 
will be heated up using the blackfold approach. At the same time we highlight that, beyond 
the usual 1 /2 BPS solution, there is a stable branch of giant gravitons that has not received 
much attention in the literature. Some unnoticed properties of this branch will be discussed 
as well. 

2.1 Setup and action 

We consider ten-dimensional type IIB string theory on AdSs x S 5 with radii L and five-form 



flux 



F (5) = 4dfi (5) /L 



(2.1) 



where dil^ is the unit volume form on the 5 5 . For the S 5 we take the parameterization 



dttf 5) = L 2 dC 2 + cos 2 C#i + sinC 2 dft 2 . 



(3) 



(2.2) 
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where d0 2 3 ^ denotes the line element on S 3 (with coordinates 4>2,(j>$,9). The giant graviton is 
obtained by considering a (rotating) probe D3-brane in this space, that wraps an S 3 inside the 
S 5 . Denoting the world volume coordinates of the D3-brane as {<t° = r, a 1 }, its embedding 
into the background is taken to be 

t = t, 4>i = fir, 4>2 = o~i, 4*3 = &2, = a 3 , C = const. , (2.3) 

while the D3-brane sits at the origin of the AdSs space. The size of the giant graviton is 
thus r = Lsin£ and this configuration rotates with angular velocity Cl on the <S 5 , satisfying 
the geometric bound (L 2 — r 2 )f2 2 < 1. The resulting induced metric on the D3-brane world 
volume is 

-f ab da a da b = -k 2 dr 2 + r 2 dQ 2 {3) , (2.4) 

where a = r, 1, 2, 3 runs over the world volume directions and 

k = = ^1-£1 2 {L 2 -r 2 )' (2.5) 

is the norm of the rotational Killing vector satisfying k < 1. 

With this setup, the giant graviton is found by solving the equations of motion (EOMs) 
of the D3-brane DBI action Jdbi m this background. Defining the corresponding Lagrangian 
via Ibbi = / drL-DBi, we have 



LdBI = £ = ~Td3 / (V-J ~ Arai^ffa) > ( 2 - 6 ) 

Js 3 Js 3 

where 7 is the determinant of the induced metric (2.4) on the D3-brane and A TaiU2a3 the 
pullback of the four-form gauge potential onto the world volume. Using the embedding above 
this gives 

Ldbi = -T D3 ft ( 3)r 3 (k - rO) . (2.7) 
The angular momentum and Hamiltonian are then computed as 

( ,s, 

We finally note that the overall factor in all these expressions involves Td3^(3) = N/L 4 where 
N is the background flux. 

2.2 Solution branches and stability 

Varying the Lagrangian (2.7) with respect to r we obtain the EOM 



3 - 3L 2 Q 2 + Arn (rn - ^l-Q 2 (L 2 -r 2 ) ) =0 . (2.9) 



This equation has two branches of solutions 
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= ^ + = V9L2 3 _ 8r2 ^ (0<r<L), (2.10) 



3 The limit r — of these solutions describes the point-particle limit of the giant graviton, where one should 
be careful in taking the limit r — > such as to obtain sensible conserved charges [15]. 
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which we call the lower and upper branch respectively. Note that for the upper branch we 
have that 1 < < 3. These two solution branches are depicted in the left plot of Fig. 1. 
It is interesting to note that a maximal size giant graviton (r = L) exists in both branches 
with either Cl- = L _1 or £l + = 3L _1 . Moreover it is also worth noticing that both branches 
connect to the point-particle case in the limit r — > 0. 



(l J 




Figure 1: = 0,L (left plot) and J = J/N (right plot) versus f = r/L for the two solution 
branches of extremal giant gravitons. The lower (— ) branch is blue and the upper (+) branch 
is red. 

To elucidate these branches and connect to a more physical parameterization we use (2.8) 
to compute the on-shell angular momentum J(f) and energy E(f) for each of the two branches, 
where we introduce the dimensionless ratio r = r/L, yielding the angular momentum 

J_=Nf 2 , J + = Nf 2 (3-2f 2 ) , (2.11) 

and energy 

E- = ^-f 2 , E + = —r 2 \/9-8r 2 , (2.12) 

L L 

on each of the two branches. For clarity, we have depicted these results in the right plot of 
Fig. 1 and the left plot of Fig. 2 respectively. Here and in the following we have also rescaled 
J = J/N and E = (L/N)E. 

One observes that for each value in the range < J < 9/8, there are two possible solutions, 
with different values of r. Comparing the two corresponding values of the energy for each of 
these two values of r (given J), one finds that the one with highest f minimizes the energy. To 
see this more clearly, we exhibit E versus J in the right plot of Fig. 2. Hence we expect that the 
stable branch of solutions consists of the entire lower branch (for < J < 1 and < f < 1) 
together with the part of the upper branch that has 1 < J < 9/8 and Vo/2 < r < 1. 
Conversely, the part of the upper branch spanned byO<J<9/8 and < f < '/2 will be 
for given J a local maximum of the energy. 
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More properly, this result on the dynamical stability can be derived by computing the 
off-shell Hamiltonian from (2.8) 



Varying this with respect to r for constant J gives, as expected, the extrema SI = Cl± found 
before. To see which part of the branches are stable we vary H once more with respect to f at 
constant J, and demand positivity, so that we are at a minimum. The result is that the lower 
branch SI = Jl_ is stable for all values of r (0 < r < 1) and the upper branch = Q + is stable 
for v / 3"/2 < r < 1. This is in accord with the arguments of the previous paragraph (see also 
App. A where the same conclusion is obtained from a more detailed stability analysis that 
includes time derivatives of the radial coordinate). Finally, we note that the point r = \f?> jl 
where the upper branch becomes unstable can also be seen as a turning point in a plot of f2 
as a function of J (see Fig. 3). 

The main motivation of our review above and the various plots that are presented is that 
they will be instructive to illustrate the new features that appear when constructing and 
analyzing the thermal giant gravitons in Sees. 4 and 5. 



Having established which solutions are stable, we now turn to their physical relevance. 
First we note that they are distinguished by the angular momentum J. In terms of r they 
coexist when \/3"/2 < r < 1, but from (2.12) we easily see that the energy on the lower 
branch is lower than that of the upper branch, for given r in that range, except when f = 1 
where they have the same energy. The lower branch is the usual 1/2 BPS branch extensively 
considered in the literature, and from (2.11), (2.12) we immediately see the BPS condition 
E- = J-/L. The other stable solution which is part of the upper branch was noted in [15] 
(see in particular Figs. 1 and 2 of that reference), but has otherwise largely been ignored. 
First of all this branch is not connected to the point particle case as a stable configuration 
since local stability requires \[?> jl < r < 1. Furthermore, while it is a local minimum of the 
energy it is not a global one, so it is a metastable configuration and has E+ > J+/L where 
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Figure 3: f2 versus J for the two solution branches of extremal giant gravitons 

the bound is saturated for r = 1. This thus raises the question whether this configuration 
indeed preserves 1/2 of the supersymmetries. By repeating the steps of Sec. 3 of [15], we have 
verified that this is indeed the case. The main point here is that the O-dependent terms in 
this computation vanish at r = 1. So we see that at f = 1 we can have either fi_ = 1/L 
or = 3/L, both satisfying the same BPS bound and both being supersymmetric. In 
particular, we cannot distinguish these configurations according to their energy. 4 In the first 
case the center of mass is rotating at the speed of light while in the second the center of mass 
is rotating at a superluminal velocity. However this should not be an argument for discarding 
the latter solution since the center of mass being a geometrical construction can be moving 
with superluminal velocities as long as every point on the brane is subluminal. The existence 
of these two BPS configurations at f = 1, arising from two distinct solution branches raises 
the question of what the dual CFT interpretation is of the one connected to the non-BPS 
branch, and we briefly comment on this in the next subsection. 

The primary purpose of this paper though is to use the blackfold approach to find the 
thermal versions of the extremal configurations reviewed in this section. Concretely, we will 
switch on a temperature and examine what happens with the entire solution branch depicted 
in Figs. 1, 2 and 3 irrespective of their stability properties and subsequently identify the stable 
branches. As it will be discussed in more detail below, in the blackfold approach one obtains 
solutions in the supergravity (closed string) regime. In this sense the entire solution branch 
(stable or not) is a valid approximate analytical solution of the type IIB supergravity EOMs, 
at least in the pertubative regime where the approach is valid. Moreover, just as a subset of 
the branches in the extremal case is locally stable, we will likewise see that a corresponding 
subset is dynamically stable when the system is considered at non-zero temperature. 

2.3 CFT dual and correlation functions 

The CFT dual operator of a single point-like graviton is a chiral primary of the form 

C = TrZ J , (2.14) 
4 They also satisfy the zero temperature limit of a general Smarr relation that is derived in App. B. 
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with J the angular momentum on the S 5 and Z a complex scalar field. Standard computations 
have shown that their two- and three point functions match exactly on both gauge and string 
theory sides provided J is small. 

If J S> N/y/X" the correct description is in terms of a giant graviton. The dual gauge 
theory operator O gg of the giant graviton is no longer given by (2.14) and arguments based 
on symmetry (which only really apply close to r = L) imply that it must be replaced by a 
Schur polynomial operator of the form [17] (see also [18]) 

O sg ~ XR{Z) = ± J2 XrA^Z^ ...ZZ {J) , (2-15) 



res. 



where Z is a complex matrix, R n denotes an irreducible representation of U(N) described in 
terms of a Young tableau with J boxes. 

As explained above, there is another (upper) branch of giant gravitons which is 1/2 BPS at 
r = L in the large J limit with the same quantum numbers as the lower branch. We speculate 
that there exists another 1/2 BPS Schur polynomial operator in the CFT at J = N that is 
distinct from the Schur polynomial relevant to the usual (lower) BPS branch and which is 
dual to the upper branch of giant gravitons at r = L. We present indications of this below. 

Two-point correlation functions 

As an explicit check of the statement above, we now compute the two-point function for the 
CFT operator dual to the r = L point on the upper branch, showing that it has the same 
properties as the r = L solution of the lower branch. It is easiest to do the computation 
simultaneously for both branches. Our method is based on the general prescription, reviewed 
in [19], for computing two-point correlation functions for massive (or light) particles moving 
in a background spacetime. 

The giant graviton is a brane, not a particle, however as seen from the AdS5 part it is 
a point-particle with a certain mass [20]. This can be seen by introducing motion in the 
AdSs part, i.e. introducing the dependence x^(t), /j, = 0...4 on the coordinates of AdSs 
with metric G^. Following [20] one can then show that the DBI action can equivalently be 
written as 

1 f (Gn V ±^x u Q 2 (L 2 -r 2 ) » 2 \ , 
/dbi =2 dT { e + ~^~ e " m 6 + m rQ ) ' (2 " 16) 

where we have defined m = Nr 3 /L 4 and e is an einbein which acts as a Lagrangian multiplier. 
Using (2.8) we can eliminate in favor of J and arrive at the action 
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where we have defined 



J 2 — L 2 m 2 

M =\l L 2_ r 2 ■ ( 2 - 18 ) 

However, to arrive at the interpretation that from the AdSs perspective the giant graviton 
is a massive point particle moving along a timeline geodesic, one should take into account 
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that J must be conserved along any path. Hence, one should consider the Routhian 1Z which 
is obtained by doing a Legendre transformation in the cyclic coordinates. In this case it 
coincides with the Hamiltonian, and hence we find 

K = H = nj-L = -± ( G^x v _ &E 2^ (2ig) 

where E is the on-shell energy (2.12) for each of the branches. So we find that from the AdSs 
perspective the giant graviton is a point particle with mass E. Following [19], wee can now 
compute the two-point function using the Routhian 



x 



-2E± 



G(0, e; x, e) = e~ K ~ ( g ) > ( 2 - 20 ) 

showing for both branches equality of the anomalous dimension and the energy. We thus con- 
clude that the anomalous dimension of the operator is equal to the energy for both branches, 
thus giving strong indication of being in both CELS6S £1 Schur polynomial at the r = L point. 

It is important to note that the correct result is reproduced here using the Routhian, and 
not the action, as was also advocated recently in [21]. Indeed, evaluating the quantity M in 
(2.18) for each of the solution branches found in subsection 2.2 one finds 5 

M_ = ^-f 2 = £L , M + = — f 2 V9 -if 2 ' / E + , (2.21) 
as compared to the energies given in (2.12). 

Three-point correlation functions 

To gain further insight into the nature of the new state at r = L one may consider the 
three-point correlation function between one point particle and two giant gravitons. For 
the lower branch this analysis was performed in [22]. The procedure consists in analyzing 
the supergravity modes which describe fluctuations in the Euclidean D-brane action of the 
metric and 4-form potential, which are dual to chiral primary operators with R-charges in the 
N = 4 SYM theory. The resulting three-point function structure constant for the maximal 
size 1/2-BPS giant graviton was found to be zero in agreement with the gauge theory side. 
Following the same steps for the upper branch r = L state gives zero as well, since one can 
check that in that case the result is independent of f2. This provides further confirmation that 
the gauge theory description of the upper branch r = L state is a Schur polynomial. It would 
be very interesting to calculate this three-point function more generally for the entire (non- 
BPS) upper branch, but this is beyond the scope of the present paper. A naive application 
of the ideas mentioned above does not give sensible results, so perhaps one should use the 
Routhian rather than the action and/or introduce an appropriate cutoff to regularize the 
divergent integrals. 



5 In Ref. [22] the action was used to compute the two-point function, but since this computation was for the 
lower (1/2-BPS) branch, for which the terms conspire to give M_ = E-, this still gives the correct result. 



10 



3 Blackfold method for branes in flux backgrounds 



In this section we briefly discuss how to extend the blackfold method for branes embedded in 
a background with fluxes. Moreover, we present a slight generalization of the previously ob- 
tained results for stationary blackfolds to blackfolds that are "boosted stationary" blackfolds, 
with the boost along a Killing direction of the background. These ingredients are relevant to 
the application of the blackfold method in Sees. 4-6 to construct and analyze localized thermal 
giant gravitons based on black D3-branes wrapping an S 3 in the AdSs x S 5 background. 

Our conventions follow [2, 6] (see also [23]). In particular, space-time coordinates are 
denoted by x^, brane world volume coordinates by a a . Furthermore, the brane embedding in 
the background space-time with metric g^ u is denoted by X p {a), so that ^ a b = g^daX^dbX 1 ' 
is the induced metric on the brane world volume. We will also need the orthogonal projector 
U* v defined by g pv = h» v + where the tangential projector is bT = ~i ab d a X»d h X u . 

3.1 Blackfold EOMs 

Consider a p-brane embedded in a D-dimensional background. In the (probe) blackfold ap- 
proximation we regard this brane as infinitely thin and thus having a localized stress-tensor 

f^(x) = f d p+1 a^Y 6{D){x Z^ {a)) T^(a) , (3.1) 
J V=5 

where the integral is over the world volume W p +i of the brane. In the following we furthermore 

consider the background to have a (p + 2)-form flux F^ p+2 ) = (L4( p+1 ) under which the p- brane 

is electrically charged with charge Q p . The charged p-brane therefore has the (p + l)-form 

current 

J {p+ i)(x) = J de +l o^ 5{D){x ^ {a)) J (p+l) {o) , (3.2) 

with 

J(p+i) = QpU( p+ i) , W(p+i) = da A • • • A da p . (3.3) 

The EOMs for the charged brane probe in a background flux are 

V ^ = frhv FUpi ''' Pp+lJpv -- p * +i ' (3 - 4) 



v ^p Vl ... Vp = ^ ^3 ^ 

where the righthand side of (3.4) includes the generalized Lorentz force. By projecting these 
equations on directions parallel and perpendicular to the world volume of the brane, we 
get a set of intrinsic and extrinsic equations respectively. The intrinsic equations include 
D a j aa i--- a P = Q ; which correspond to charge conservation as well as 6 

D a T ab = , (3.6) 

6 In the case of a submanifold with boundaries (with unit normal vector h a ) these must be supplemented 
by the boundary conditions T ab hb\aw p+1 ~ an d J a n a \dw p+1 = 0. 
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under the assumption that the background gauge-potential is constant along the world- volume 
of the brane, so that there are no forces induced on the world volume. This will be the case 
in our applications. The extrinsic equation for the brane embedding is (see e.g. [24, 2, 6]) 

TabR ^ = ^jlV 1 "^ 1 ^...^ , (3.7) 

where 

K ab =^x {dad b X x + T x up d a X»d b Xp} (3.8) 

is the second fundamental tensor. 

We note here that the EOMs for extremal F-strings, D/M-branes in backgrounds with 
appropriate non-zero flux and constant dilation field are of the form (3.7) with appropriate 
T ab . In particular, the D-brane equation of motion coming from the DBI action (including the 
WZ term), is of the form (3.7) with T ab = —T£) P ^ ab in the case of zero world volume gauge 
field. In the blackfold approach this energy momentum tensor is replaced by the effective 
stress tensor of the corresponding non-extremal brane solution of supergravity. In the leading 
order approximation this is that of a perfect fluid in local thermodynamic equilibrium, with 
equation of state determined by the brane solution. Solution of the intrinsic equation (3.6) 
requires in that case that the fluid velocity u a = k a /\k\ is aligned with a timelike world volume 
Killing field and following [2] we assume that this Killing field k a = d a X^(a)k p can be pushed 
forward to a corresponding timelike background Killing field k^. This will be the context in 
which we will derive conserved quantities below. 

Stationary and quasi-stationary solutions 

From the above, it follows that in general the normal to the spacelike hypersurface B p of 
the blackfold world volume W p +i need not be parallel to the generator of asymptotic time 
translations, i.e. one can have 

d T = < + b X , (3.9) 

where £ corresponds to the canonically normalized generator of time translation and x that 
of a spatial U(l) isometry of the background. However, with regard to stationarity there is 
an important further distinction to be made, depending on whether: 7 

i) x is a l so a worldvolume Killing vector 

ii) x is perpendicular to the world volume (and hence not a world volume Killing vector) 

In case i) the blackfold world volume does not break the isometries £, x of the background 
and in particular the conserved quantities associated to these are also conserved for the entire 
solution consisting of background with the blackfold in it. 8 The resulting solutions are sta- 
tionary blackfolds. In case ii), which is the one relevant for the present paper, the blackfold 
world volume only preserves a particular combination of the isometries £, x- As a result the 

7 For simplicity we restrict to these two possibilities, but we note that when \ is not a world volume Killing 
vector, it could also have some components along the world volume. E.g. \ could be a linear combination of 
a woldvolume Killing vector and a perpendicular component along an isometry of the background. 

8 This is for example the case for the configurations considered in Ref. [25] using the blackfold method. 
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conserved quantities associated to these are of a different nature, namely they refer to quanti- 
ties for the probe blackfold in the background space time but are not seperately conserved for 
the background including the blackfold. Only an appropriate linear combination is conserved, 
according to (3.9). In particular, the conserved quantity generated by £ should be thought of 
as the total energy E (so not the rest mass of the object) and the quantity generated by x as 
the transverse momentum J corresponding to the boost. So we see that in this case the black- 
fold is transversely boosted along a Killing isometry of the background, and hence it should 
be viewed as a "boosted stationary" solution. We will refer to this below as quasi-stationary, 
since we still have that, seen from the world volume, the blackfold is independent of time. 
This is in fact precisely what happens for the localized giant graviton, since in that case x 
lies in a direction perpendicular to the world volume. It is important to note that since the 
quasi-stationary blackfold is not accelerating it does not emit radiation and one can thus go 
beyond the probe approximation and perform a matched asymptotic expansion for the full 
system of the background with the brane. 

3.2 Conserved charges 

We now write down the expressions for the conserved charges corresponding to the asymptotic 
generators £ and x f° r these quasi-stationary blackfolds in flux backgrounds. Note that the 
results below can also be used for stationary blackfolds in flux backgrounds. 

For any Killing vector field (KVF) k of the background we have by definition that the Lie 
derivative along k of the (p + 2)-form F = F^ p+2 ) is zero C k F = 0. Since dF = we find that 
= dF = ikdF + d{ikF) = d{ikF) where it means the contraction with the KVF k. Picking 
a gauge in which C k A = we see indeed that = C k A = i k F + d{i k A) thus i k F = —d(i k A). 
Thus, in this gauge the (p + l)-form i k F has the p-form potential i k A. Using this with the 
EOM (3.4) and current conservation (3.5) we see that 

p\ V^TTK) = ^k»F upi ... Pp+1 J^-^ = -V^{i k A) P2 ... Pp+l] J^-^ 

= -V Pl (i k A) P2 -.. Pp+1 J p ^ +1 = -V fl (A vpi ... Pp J^-^) . (3.10) 
Thus we obtain the conserved current 

# = (?T + ^Pv-pJ^-^ku , (3.11) 

and a conserved charge 

Qk= dx D ' 1 ^/g spSLCe j^n pj , (3.12) 
Jt, 

where g S p aC e only involves the spatial components of the background metric, defined by the 
slice S of constant x° = t and is the unit normal of S. Inserting the conserved current 
(3.11) and using the form of the stress tensor (3.1) and current (3.2) we can do the (^-function 
integrals to reduce to integrals over B p . On the world- volume we choose the static gauge in 
which X° = t = do. Restricting to static backgrounds, which is sufficient for the applications 
in this paper, we have \J—g = y 7 — 5oo -y/Sspace and since we assume that the pullback £ a 
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is hypersurface orthogonal, we split \/ — 7 = - v /— ^y TT dV p . Then, integrating out the delta 
function we obtain for the conserved charge 

Qk= [ dV^l 1 [T^ + VH n^K\ x , =x , . (3.13) 

Here we have defined 

V^ = ^ A\ 1 ,.. tlp _ 1 J^--^ , (3.14) 
and the transverse Lorentz contraction factor 7j_ is given by 

7± = - , (3-15) 
Po 

where i? = yJ-Qtt and p = \I~1tt- 

Now, we use the result (3.13) to write down conserved charges corresponding to the back- 
ground Killing vectors £ and % 

E = [ dV {v)^l l + V H > J = ~ [ ^(rtTl 1 + VH n M x, , (3-16) 

which are the energy E and momentum J of the quasi-stationary blackfold moving with 
constant velocity in the background along an isometric direction. The interpretation of these 
conserved charges in this case is most easily seen by analogy with a probe particle moving 
in a time-independent background along a Killing direction. In that case, it follows from 
standard analytical mechanics that the energy E and momentum P of the probe particle are 
conserved as long as the object moves with constant velocity. Likewise for the blackfold, as 
long as we are working in the leading order probe approximation where the laws of physics 
do not involve the internal degrees of freedom, the conservation of E and J relies entirely 
on properties of the background and not those of the brane. In particular, these quantities 
are conserved for the quasi-stationary blackfold probe just as they are for the DBI D-brane 
probe. The expressions (3.16) will thus provide us the conserved quantities for the thermal 
giant graviton constructed in Sec. 4, where we will also see that they reduce to the correct 
DBI quantities in the extremal limit. 

3.3 Action and thermodynamics 

Finally, we discuss here the action that describes the quasi-stationary blackfolds introduced 
above. Again, we emphasize that the considerations of this section can also be applied to the 
stationary case. 

When the intrinsic EOMs are solved as explained in Sec. 3.1 the remaining extrinsic EOMs 
(3.7) can be shown to follow from the Lorentzian action 

1=1 {u {p+ i)P + Q P nA {p+1) }} , (3.17) 

where P is the local pressure of the blackfold and "PLA^+i)] the pull-back of the background 
gauge potential to the worldvolume. This natural generalization of the blackfold action to 
include background fluxes, easily follows form the derivations presented in Refs. [2, 6]. For 
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a given set of global parameters (T, f2, Q p ) the (quasi)-stationary solution to the blackfold 
equations of motion is exactly the one that extremizes the action 51 = 0. Note also that for 
extremal D-branes this action reduces to the DBI action plus WZ term. 

Following the arguments in Refs. [2, 6], we show in App. B that the Euclidean action Ie 
obtained by Wick rotating (3.17), is again equivalent to the thermodynamic action 

= F = E -07 -TS . (3.18) 

Here the total entropy S of the blackfold is obtained in the usual manner by integrating the 
temporal part of the entropy current su^ over B p , so that 

S = I Tl 1 ^^ , ( 3 - 19 ) 

where the integrand is multiplied by the appropriate Lorentz factor 7j_ defined in (3.15). 
Consequently we find that (at fixed Q p ), the extrema of the action obey the first law of 
thermodynamics dE = fid J + TdS. We finally note that a corresponding Smarr relation for 
the thermodynamic quantities is derived in App. B as well. 



4 Construction of finite temperature giant graviton on S 5 

In this section we find a thermal version of the giant graviton consisting of a D3-brane wrapped 
on a 3-sphere moving on the 5-sphere of AdSs x S 5 (reviewed for the extremal case in Sec. 2). 
This is done using the blackfold approach in the probe approximation (as reviewed in Sec. 3.1). 

4.1 Equation of motion for thermal D3-brane giant graviton 

In the leading order blackfold approximation to black branes we use the general extrinsic 
equation (3.7) with the stress-tensor T a ^ being that of a black brane. In particular, in the 
present case we want to consider black D3-branes in the AdSs x S 5 background. Specifically, 
the stress-tensor of black D3-branes corresponds to that of a four-dimensional fluid tensor of 
the form T a b = (e + P)u a Ub + P^fab with u a being the four-velocity and energy and pressure 
given as 

vr 2 

e = Ts-P , P = -g (1 + 4sinh 2 a) , Ts = AG , Q = — T^ri (4.1) 
where the local temperature T and entropy density s for the black D3-brane are 

T = — — , s = 27r 3 r^ 3 r^ cosh a (4.2) 

7rro cosh a 

The parameters of the black D3-brane stress tensor and thermodynamics are thus ro and a. 
The black D3-brane furthermore has the 4-form charge current 

^(4) = Qdrda da 2 da , Q = 4(/ sinh a cosh a = Nd^Td3 (4-3) 

where Q is the charge density and A^3 is the number of coincident D3-branes and Td3 = 
l/((2ir) 3 g s lg) is the D3-brane tension. Note that we have the relation S7( 3 )Td3 = N/L 4 where 
and L are the magnitude of the flux and the radius in AdSs x S 5 , respectively. 
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The giant graviton that we wish to thermalize is that of a D3-brane wrapped on a 3-sphere 
moving on the 5-sphere of AdSs x S 5 . This giant graviton solution, as studied using the DBI 
action for the extremal D3-brane, was reviewed above in Sec. 2. In particular, the metric on 
the 5-sphere is given by (2.2) and the 5-form flux in (2.1). We now use the properties (4.1)- 
(4.3) of black D3-branes to study the case where the D3-brane is put at finite temperature, 
hence we are studying what can be called a thermal giant graviton in the form of a black 
D3-brane wrapped on a 3-sphere moving on the 5-sphere of AdSs x 5" 5 . We take the same 
ansatz for the embedding (2.3) as used for the extremal D3-brane. Using all this information 
(see App. C for details of the derivation) 9 with Eq. (3.7) we arrive after some algebra at the 
following EOM 

Q 2 r 2 (1 - fti(a)) + 3k 2 + 4MlrK 2 (a) = , (4.4) 



where k = \k\ = y/l — Q 2 (L 2 — r 2 ) and we introduced the following quantities 

Ts 4 Q 4sinhacosha 

^ l(a) "^ = -l + 4sinh 2 a n ^^P = - l+4sinh 2 a ' (4 " 5) 

We notice that by taking the limit a — > oo we obtain the EOM (2.9) for the extremal case. 

Since the D3-brane is moving on the 5-sphere the local temperature has a redshift as 
compared to the global temperature T of the background space-time that we are probing as 
T = T/k. Thus, we are imagining finding the thermal giant graviton solution for a given value 
of T. In addition, the above EOM (4.4) should be supplemented by the charge quantization 
condition (4.3) which from the above becomes 

^ 3T4 = ^^iuha k4 _ 
tt cosh a 

Our goal is to study the thermal giant graviton for a given value of A^D3, T and J, or, 
alternatively, N^, T and r. From the latter choice we see that a = a(Nj)3,T,r). 

Regarding the validity of our black D3-brane probe in AdSs x 5" 5 we see that we need 
Ad3 3> 1 to have a valid SUGRA solution of the probe but that at the same time we also 
need <C N since the back reaction of the Ad3 D3-branes of the probe should be negligible 
in comparison to the back reaction of the N D3-branes from which the AdSs x background 
originates. 



4.2 Solution space 

In this section we examine the structure of the solution to the EOM (4.4) and the charge 
quantization constraint (4.6). 

The aim in the following is to study the thermal giant graviton solution given a temperature 
T, the number of D3-branes A^D3 for the brane probe as well as the radius r of the 3-sphere 
that the D3- brane probe is wrapped on. 10 We find the following two branches of solutions 

9 Note that the worldvolume velocity field is given by u T = 1/k so on the worldvolume the D3-brane is 
static, but the push-forward of this vector to the background gives the vector field it M 9 M ~ dt + £ld$, so this is 
a quasi-stationary blackfold as explained in Sec. 3.1. 

10 One can alternatively use (5.1) to write r — r(A r D3,T, J) and consider Nr>3, T and J given, which will be 
further examined below. 
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Ofc(r) to the EOMs (4.4) 



n± to = , — ; , (4.7) 

'9L 2 -8(1 + A±(a))r 2 



which we refer to as the lower (— ) and upper (+) branch respectively. Here we have defined 11 

A±(a) = -- (3^i (a) + 8ft 2 (a) 2 ± 472 2 (a)V©(a) 1 ) + ^ , 2?(a) = 3fti(a) + 47e 2 (a) 2 - 3 , 

(4.8) 

where 72-1 2(0:) are given in (4.5). These solutions are supplemented with the charge quanti- 
zation constraint (4.6). In order to determine a± as a function of T, Nds and r one should 
substitute (4.7)-(4.8) into (4.6). 

Taking the extremal limit a — > 00 we see that 7Z± — > and 72-2 — >• — 1. Hence the above 
solutions (4.7)-(4.8) reduce to the extremal values given in (2.10), found in Sec. 2.2 from 
the DBI analysis of the extremal giant graviton solution. 

Considering Eq. (4.6) we have that sinha/ cosh 3 a is bounded from above with maximal 
value 2y / 3 1 /9 corresponding to the value a for which cosh 2 a = 3/2. Thus, since we choose to 
be on the branch connected to the extremal D3-brane we always have a > a. Setting Q = 
and using this bound in Eq. (4.6) gives rise to the maximal temperature for coincident 
black D3-branes given by T s t a tic = (^f"^ 3- ) 1 ^ 4 - However in the case under investigation - 
where > and depends on a - one obtains a stronger bound on a from the requirement 
that T>(a) in (4.8) should be always non negative. Indeed, T>(a) = for a = a = cosh" 1 (3/2) 
and we see that a > a. Thus for the finite-temperature giant graviton we have the bound 
a > a. This means that from (4.6) we have the bound 

k > f , (4.9) 

where we have introduced the rescaled temperature 

*=:£-r-= pSf' <«<» 

Moreover we have the geometric upper bound r < L on the size of the giant graviton, which 
means that k < 1. It thus follows that for the giant graviton on 5 5 the value of k must lie in 
the range 

f<k<l (4.11) 

and as a corrollary we see that the temperature is bounded, T < 1. In particular, when the 
upper bound on the temperature T max is reached the solution space collapses to a point. This 
upper bound on the temperature is more restrictive than for a static black D3-brane where 
the maximal temperature is T s t a tic given above. 



Note that — > —CI, a — > —a is also a solution. Here we only consider £1 > (=> a > 0). 
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Parameterization of thermal giant graviton solution 

We now describe a very useful analytic parameterization of the solution, which will be em- 
ployed in this and the next section to analyze the solution in more detail. For given T and 
iVD3 the solution can be parameterized by the value of k, as follows. Introducing 

1 (4-12) 



cosh 2 a 

one finds that Eq. (4.6) can be rewritten as the cubic equation 



b 3 -^ 2 + ^ S m 2 5 = , (4.13) 



where we have defined 



Eq. (4.13) is of the same type as encountered for the thermal Blon [8, 9] and the solution 
connected to the extremal giant graviton is 




2 sin 5 



^T^) = -^ j — . (4.15) 

3 a/3 cos I — sin 



We thus have an explicit functional expression for a(T, k) = arccosh^>(T, k) -1 / 2 and substi- 
tuting this in (4.8) we then obtain A±(T, k). With those expression in hand, one can now 
obtain r as a function of (T, k), using k = \/l — fl 2 (L 2 — r 2 ) and the solution (4.7). This 
yields 

3k 

r±(T,k) = — f —L , (4.16) 

V8k 2 (l + A ± (a)) + l-8A±(a) 

along with 



where we recall the range (4.10) for k. 

Fig. 4 depicts the resulting solution branches for various values of T in a (f = r/L, $7 = VtL) 
diagram, and we recall for comparison that the corresponding extremal solution is plotted in 
Fig. 1. We note the following new and interesting features. The lower and upper branch meet 
in the point where k = T (i.e. a = a) saturating the lower bound in (4.9). At this point we 
find A± = — \ giving Sl± = Q, = 3/\/9L 2 — 4f 2 . Inserting this into the expression for k in 
(2.5) then gives 

r = , 3L , • (4.18) 
V4 + 5T - 2 

At the upper bound k = 1, we have r = L and S7_(T, 1) = 1/L and fi + (T, 1) < 3/L. 
Furthermore, we observe that the values of r are restricted to < r m i n (T) < r < L, and 
fm.in{T) approaches L as the maximum temperature is approached. The minimal size thermal 
giant graviton r m j n lies on the lower branch, which curves back to meet the upper branch in 
the point f. Furthermore, for each r value in this range there are two possible values for O, 
lying in between the two corresponding values of the extremal solution. 
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Figure 4: O = versus f = r/L for the two solution branches of thermal giant gravitons for 
various values of T. The lower (— ) branch is blue and the upper (+) branch is red. Shown 
are the values T = 0.1, 0.2, 0.3 and 0.4 with temperature increasing as the curves move to 
the right and become fainter in color. 

That the minimal size of the giant graviton is greater than zero is an important consequence 
of the finite-temperature physics of the giant graviton. For the extremal giant graviton the 
two branches meet in the singular solution r = which in turn corresponds to the graviton 
particle with same angular momentum. What we see at finite temperature is that: a) there 
is a minimal possible size r m [ Q of the giant graviton and b) unlike in the extremal case, it is 
possible to move in the solution space from one branch to another since the meeting point of 
the two branches at f is not a singular solution. Note that the fact that the thermal giant 
graviton attains a minimal possible size has an analogue in the thermal Blon as well as in the 
thermal Wilson line cases studied in [8, 9, 10]. 

In the next section we will use the parameterization above to make a detailed analysis of 
the thermodynamic quantities of the two thermal giant graviton branches and examine their 
stability. In particular, the explicit parametric solution enables to easily plot any combina- 
tion of thermodynamic quantities and study their behavior. However, the expressions are 
too complicated in general to invert but there are three interesting situations which will be 
considered in more detail. 

• The low-temperature regime: We can expand around the extremal case by doing 
an expansion for small T. 

• The maximal size regime: We can expand around k = 1 (i.e. r = L). This also 
includes the maximal temperature regime T — y T max since in that limit the physics is 
captured by the large giant graviton limit e = (L — r)/L -C 1. 

• The minimal charge parameter limit: We can expand around k = T, i.e. the point 
a = a where the two branches meet at the value f. 

The first regime above, which is physically the most relevant one, will be investigated 
further in Sec. 5, while the other two limits are considered in App. D. Since our main focus 
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in the text will be the low-temperature or near-extremal regime, we now present the solution 
described above in this limit. 



Low-temperature solution 



The convenient feature of the parameterization above, using cj)(a) defined in (4.12), and 5 
defined in (4.14), is that the extremal limit is obtained for <p ~ > and 5 — > 0. Expanding 
around this, we can work out the form of the angular velocity f2 and charge parameter <f> as 
a function of the thermal giant graviton size r in the low temperature limit. 

In order to do a low temperature expansion we demand that T< 1 and sin 5 <Cl. As a 
consequence, we need that k 3> T in order for the expansion to be valid, which in turn implies 
that r 3> TL. We then have to leading order in T, 



S 




3V3 



6 . 



(4.19) 



where we used (4.14), (4.15). Performing a small 5 expansion of A± gives then to leading 
order A + ~ — </>/4 and A_ ~ —1 (here and below, expressions are up to corrections of order 
T 8 ). This can now be substituted into the expressions (4.16) for r±, yielding 



kL 



kL 



Vl + 8k 2 




1 



l + 8k 2 



(4.20) 



These expressions can be inverted in order to write k as a function of r. We find 

rCl + (r) t 4v / ET p 2 



+ 



rQ + (r) r 2 



(4.21) 



where Cl± are the extremal values (2.10) and where we have defined p 2 = r 2 — L 2 , while the 
size r now parameterizes the solution and the ± subscript has thus been moved from r to k. 
Substituting (4.21) into the expression (4.17) for O we then obtain 



along with 



Q_(r) 



<f>-( r ) 



1 
L 



n+(r) ~ n+(r) 



1 



4V5 



3r 2 fi2 (r) 



^ , Mr) * 12V5 



rp4 



(4.22) 



(4.23) 



27 r 4 ' ' r 4 ^ 4 (r) ' 

Note that the angular velocity on the upper branch is much more sensitive to low temperatures 
than the lower branch. 



4.3 Validity of probe approximation 

For the probe approximation to be valid for our SUGRA black D3-brane probe we must 
require the transverse length scale r s of the probe to satisfy the following conditions 

r s < Hnt , r s < r ext , r s < L , (4.24) 
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where r- m x and r ex t are the length scales associated with the intrinsic and extrinsic curvature 
of the embedding of the brane, respectively, and L is the length scale of the AdSs x S 
background. For the black D3-brane in the branch connected to the extremal solution the 
transverse length scale r s is easily seen to be given by r s ~ (y^) 1 / 4 [8]. 

We compute the Ricci scalar for the embedding metric in order to obtain the intrinsic 
length scale. This simply gives r mt = r/^/W. Instead the extrinsic length scale of the embed- 
ding is obtained as r ext = \KPrip\~ 1 where n p is the unit normal vector to the brane embedding 
and K p is the extrinsic curvature. We find 



rk 2 7k 2 + l-r 2 /L 2 ' 
rext " (l-g) tt 2 r 2 + 3k 2 ■ ( 5) 

Collecting now this information with (4.24) we see that we have two different regimes to 
consider for the validity of the probe approximation, namely whether r/L is small or not 
(note of course that < r/L < 1). If r/L is not small, then we can roughly regard r and L to 
be of the same order, hence (4.24) simply reduces to r s <C L which using L 4 = N/(Q,^T-£)3) 
one can write as 

N m « N , (4.26) 

a condition already mentioned in the beginning of this section. Instead, if r/L is small we 
should impose the condition that r s <C r which we can write as N^/N <C r 4 /L 4 . 

In addition to the validity of the probe approximation we should also require the validity 
of the SUGRA description of the black D3-brane probe. This requires A?D3 S> 1 as well as 
9sNy>3 ^> 1- The latter condition can be written as AA^D3 3> N (A being the 't Hooft coupling). 
Therefore, if we assume we are in the regime where r/L is not small we can summarize both the 
probe approximation condition as well as the conditions for validity of the SUGRA D3-brane 
description as the conditions 

1 <C iV D3 < N < XN D3 . (4.27) 

Hawking-Page temperature 

An important point is to examine how the bounds above relate to the Hawking-Page temper- 
ature, above which the AdS black hole background will become dominant over the hot AdS 
spacetime background considered in this paper. Using that the Hawking-Page temperature 
Thp ~ and the expression for T max in (4.10), we find 



T I N 

- 1 max I ly 



1/4 

» 1 , (4.28) 



where we used (4.27) in the last step. Thus in the regime where the probe blackfold ap- 
proximation is valid the maximum temperature that the solution exhibits is far above the 
Hawking-Page temperature. Consequently, this maximum temperature is not physical in the 
sense that before reaching it one should change the background to the AdS black hole. In 
particular, this means that our solution should be considered for small temperatures (much 
less than T max ) only. This is very similar to the case of thermal string probe in AdS considered 
in Ref. [10]. 
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5 Thermodynamics and stability properties 



In this section, we further investigate the thermal giant graviton solution obtained in Sec. 4. 
We will first compute the relevant conserved quantities and thermodynamic properties using 
the formulae derived in Sec. 3.2. This will also enable us to show that the same solution 
can be derived from an action and verify the first law of thermodynamics. Then the solution 
parameterization of the previous section will be used to examine the detailed behavior of the 
thermodynamics and determine which part of the solution branches are stable. Furthermore, 
we make a detailed analysis of the low temperature regime. 

5.1 Thermodynamic quantities and first law 

Using the expressions in Eqs. (3.16), (3.19) together with the perfect fluid blackfold stress 
tensor (see (4.1), (4.2)) and current (4.3), the (off-shell) energy E, angular momentum J and 
entropy S are computed to be 

E(r) = ^ ( 2 r3 o v > J ( r ) = toE(r)(L* - r 2 ) + tt (3) Qr 4 , S(r) = n {3) r 3 s(r) , (5.1) 
Y 1 — ii (iv — r z ) 

with 

fM .Zg» (l-^ 2 -r 2 )) 2 5 + 4sinh>(r)) (l - tt 2 (L 2 - r*)f 2 

e[ ) ~ 2 7r 2 T 4 cos h 4 (a(r)) ' S[ ) ~ D3 tt 2 ^ cosh 4 (a(r)) ' (5 j 

where Q( 3 ) = 27r 2 . 

We can use these quantities to compute the thermodynamic action (3.18) (see also App. B), 
which is the Gibbs free energy, 

PI E = F = E -TS -QJ = -n (3) -p (r 3 kP + r A nQ) , (5.3) 

where P is the pressure defined in (4.1). Varying with respect to r keeping T, £1 and Q 
constant, one may check explicitly that one indeed obtains the EOM (4.4), i.e. the equation 

dr dr dr 

is equivalent to the EOM (4.4). We see that this in turn means that the first law of thermo- 
dynamics is satisfied for solutions of the EOM (4.4). 

To obtain the on-shell expressions for the conserved quantities in (5.1), it is useful to define 
the rescaled quantities 

p - LE T — J C - LT ^ S ,ctC\ 

N m N ' - N m N ' S - iV D3 iV ' 

where we recall that Nd^N = Q^QL 4 . These satisfy the first law of thermodynamics in the 
form dE = TdS + CldJ, where £1 = £IL. One can also check that the Smarr relation (B.10) 
is satisfied. In the following we will also use again r = r/L to simplify the equations. Using 
(5.1) and the definitions above, one finds 

27 k 3 r3_ 

E±(T,k) = — ^^(4 + 0), J ± (T,k) = E ± ± (l-f 2 )+r 4 , (5.6) 
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S±(T,k) 



27 k 5 f 3 



± J.2 



(5.7) 



Ay/5 f 5 ' 

where 0(T, k) is determined by (4.15) and r±(T, k) for the two solution branches given in 
(4.16). It also follows that the on-shell free energy is given by 



F±(T,k) 



27 



16V5 r 4 



0(4 - 30) 



(5.8) 



These explicit results enable one to easily plot and examine any combination of the above 
quantities including f± for given value of T. 

To start, we exhibit in Fig. 5 the solution branches in a (f , J) plot for various values of T, 
and we recall that the corresponding figure for the extremal case is given in Fig. 1. We note 
the following new features. The angular momentum lies in the range J m \ n (T) < J < J m ax(T), 
where the boundary values satisfy dJ{r)/dr = 0. We will denote the corresponding r values 
by rj min and rj max respectively. Contrary to the extremal case we observe that one has a 
non-zero lower bound Jmin (T) which increases with T. On the other hand, the upper bound 
decreases with temperature. Just as for the extremal case, we observe that for each J in this 
range, there are two solutions depending on the value of r. We will shortly see that the one 
with largest r is the one that is stable. 




0.2 0.4 0.6 0.8 



Figure 5: J versus f for the two solution branches of thermal giant gravitons for the same 
values of T as in Fig. 4. 



5.2 Stability 

To address the stability we turn our attention to the on-shell free energy given in (5.8). In 
Fig. 6 we have depicted (f , F) as well as (J, F) plots for various values of T. Note that the 
corresponding plots for the energy in the extremal case were given in Fig. 2. Comparison of 
the free energies then shows that the lower branch is expected to be stable for ?"j m . n < r < L 
(with J m i n < J < 1) and the upper branch for rj max < r < 1 (with 1 < J < J m ax)- This is 
entirely in parallel with the stability properties of the extremal giant graviton (see Sec. 2.2), 
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the difference being that as a consequence of the finite temperature, a part of the lower branch 
has become unstable and there is a minimum angular momentum. Note that it follows that 
the minimum size stable thermal giant graviton is thus v j min 5 which is greater than r m i n (for 
which the solution is unstable). We also see that the point where the branches meet in f is 
always in the unstable region. On the other hand, the branches also meet in r = L, but for 
different values of Q. These special points will be considered in more detail in App. D. 



F F 




Figure 6: F versus r (left plot) and versus J (right plot) for the two solution branches of 
thermal giant gravitons at the same values of T as in Fig. 4. 

The fact that J m ; n and J max denote the onset of instability in the lower and upper branch 
respectively is further corroborated by looking at the turning points in a (J, Cl) plot, which is 
shown in Fig. 7. We see that these boundaries of stability occur precisely at the turning points 
where dJ/dQ, = 0, in accord with expectations based on the Poincare turning point method 
(see e.g. [26] and references therein). Finally, we note that these results for the stability of the 
branches are confirmed by a more detailed off-shell analysis for the three limits described in 
Sec. 4.2. The special meeting points of the two branches k = 1 and k = T, which correspond 
respectively to the maximum size thermal giant graviton and the minimal charge parameter 
solution and their stability are considered in App. D. The most interesting limit, which is the 
low temperature limit, will be considered in the next subsection. 



5.3 Low-temperature limit 

We now compute the various thermodynamic quantities for the lower and upper solution 
branch for low temperatures. Using the expansions (4.21), (4.23) in (5.6), (5.7), one finds for 
the energy, angular momentum and entropy for the lower branch in the low temperature limit 

E_(f,f)^f 2 + ^lf 4 , J_(T,f)^f 2 + ^T 4 , TS_(f,f)^^f 4 , (5.9) 
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Figure 7: ti versus J for the two solution branches of thermal giant gravitons for the same 
values of T as in Fig. 4. 

where here and in the following p 2 = 1 — r 2 . The free energy becomes 

F_(f,f)^ + ^(l-fQfS (5.10) 

and we recall that the higher order corrections in the expressions above are of order T 8 . It is 
trivial to check that the first law dE_ = (ldJ_ + TdS_ is satisfied for these quantities. 

We can also compute the on-shell free energy as a function of the lower branch angular 
momentum. To this end we invert J in (5.9) 

••-(T.J^VT+^i^. (5.11) 



Inserting this into (5.10) we then obtain 



V^ 4 



F_(T,J)=J-^r 4 + 0(T 8 ) , (5.12) 

where we note that the first correction to the extremal results is independent of the angular 
momentum. Reintroducing the dimensions from the definitions (5.5), this gives us the final 
result for the low temperature expansion of the free energy of the lower branch thermal giant 
graviton 

F(T, J)=' J L - ^N 2 3 L 3 T A + 0(T 8 ) . (5.13) 

This is one of the central results of the paper. 

Finally we can compute the ratio J/E for the lower branch. We find 

^ = L-^N 2 m (LTf + 0(T 8 ), (5.14) 

where the first term is recognized as the usual Kaluza-Klein contribution while the second 
term is due to thermal effects. 
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Repeating this procedure for the upper branch, we find 



E + (T,f) 



3f2 
O4- 



1 + 



rp4 



~f 2 (3p 2 + f 2 ) 



1 



9V5p 



TS+(f,f) ~ i-^f 4 
ni 



(5.15) 



(5.16) 



f 4 (3p 2 + f 2 )0 4 

Again it is possible to show that these quantities obey the first law of thermodynamics. 
Finally, for the free energy one finds in this case 

3f 2 / 3^' 



F+(f,r) 



(5.17) 



One may eliminate from this f in favor of J as done above for the lower branch, but the 
resulting expression involves a complicated function of the angular momentum in multiplying 
the thermal correction, so that we omit it here. 

We note that the results for J in (5.9), (5.15) explicitly show what is seen in Fig. 4, namely 
that turning on a temperature has the effect that the lower (— ) branch is pushed to the right 
while the upper (+) branch is pushed to the left. 



Stability 

We now turn our attention to stability. The method we use is described in Sec. D.3 and is 
based on an analysis of the Helmholtz free energy Fr = E — TS whose on-shell (rescaled) 
value is related to the Gibbs free energy (5.8) through the relation 



f± + nf 4 



(5.18) 



Having described the solution space for low temperatures means that we have essentially 
solved the first derivative (Fh)(i) of the off-shell free energy for T<1, r> TL. To analyze 
the stability, we thus compute the second derivative (F#V 2 ) for both branches. We find to 
leading order in T 



2f 2 



(2) 



V5(7 - 4f 2 



9r 



r-i 



and 



(Fh) 



(2) 



2Q + f 2 



4f q 



ft 2 f 4 
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40r 2 + 16r 4 ) f 4 



(5.19) 



(5.20) 



Solving for (F#)( 2 ) = determines where a solution goes from stable to unstable. Since the 
low temperature expansion is only valid for r 3> TL, we see that the entire part of the lower 
branch captured by the low temperature expansion remains stable. However, the value of r 
for which the upper branch becomes unstable is pushed up when we turn on a temperature. 



Indeed, solving {Fh , 



(2) 



0, we find that the upper branch becomes unstable at 



L 



T 4 + 0(T 8 



(5.21) 



8L 

2 " TV 3 

Here the first term is recognized as the zero temperature instability from the DBI analysis in 
Sec. 2.2. Note also that as a consistency check, the same value of r* is obtained by finding the 
maximum of J+ in the low-temperature expansion (5.16), i.e. (83 + / (dr)\ r=rt = + C(T 8 ). 
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Maximal and minimal angular momentum 



We also derive the low temperature limit expression for the maximal and minimal value of 
the angular momentum, found on the upper and lower branch respectively. The largest value 
of J is exactly attained on the upper branch where it goes from stable to unstable. So 
Jmax = J+(r*). Using this, we find 

9 Vb 



^ ^ ^ T 4 + 0(T 8 ) . (5.22) 

This expression fits nicely with the numerical data. 

The minimal value of J is attained close to r = 0. This means that an analytical expression 
for J m i n is not obtainable from the low temperature expansion (as it is only valid for f > T). 
However, we expect the following behavior for small T 

Jmin ~ . (5.23) 

It is then possible to do a fit of the numerically obtained values for J m i n . Doing this one finds 
that (3 ~ 1.89. A plot of numerical values of J m in versus T is given in Fig. 8. 




T 
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Figure 8: The minimum angular momentum J m i n of the thermal giant graviton versus the 
temperature T. 



6 Thermal giant graviton on AdSs 

In this section we analyze the case of thermal giant gravitons moving in the AdSs part of the 
AdSs x 5 5 background. We briefly review the extremal configuration of Refs. [15, 16] and 
then present its thermal generalization following the prescriptions used in Sees. 4 and 5 for 
the case of thermal giant gravitons moving on the S 5 part. 

6.1 Extremal giant graviton on AdSs 

The AdSs metric is parameterized as 

ds 2 = - (l + ^\ dt 2 + (l + ^\ dp 2 + p 2 dn 2 3) , (6.1) 
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with the four-form RR gauge field on AdS§ given by 

r 4 

At ai ...as = _ ^^(3) i (6-2) 
where the coordinates 04 parametrize the f2( 3 ) sphere in AdSs. We use the embedding 

p = r, t = r, cf>i = Qt, C = 0, a.i = Oi . (6.3) 
The intrinsic metric then takes the same form as in Eq. (2.4) now with 

k = |fc| = ^R 2 -Q 2 L 2 , (6.4) 

where the redshift factor is given by Rq = yl + "jji ■ Also in this case we find two branches 
of solutions of the EOMs 

0_ = - , 0+ = V , (0 < r < 00) . (6.5) 

-L o 

The corresponding energy and angular momentum are 

E_ = T D3 fl (3) r 2 L , J_ = T m n (3) r 2 L 2 , 

(6.6) 

E + = T m n {3) r 2 L- l {ZL 2 + 2r 2 ) , J + = T m Q {3) r 2 L^9L 2 + 8r 2 . 

In this case the angular momentum is not bounded from above. Moreover, analysis of the 
off-shell Hamiltonian shows that the lower branch is 1/2 BPS and stable while the upper 
branch is not BPS and unstable. 



6.2 Finite temperature solution 

We now examine the heated up version of the giant graviton moving on AdSs following the 
prescription employed in Sees. 4 and 5. Using the embedding (6.3) and the AdSs background 
(6.1)-(6.2) we find the EOM 

2 

^(l-^i(a)) + 3k 2 + 4k^ 2 (a) = , (6.7) 
which has the following solutions 



^9L 2 + 8(1 + A ± (a))r 2 ' 
M± = 3^2 ' t 6 - 8 ) 

with A±(a) given in (4.8). The charge quantization takes the form (4.6) with k given by 
(6.4). Note that in this case the lower bound on a is given by a with cosh 2 a = 3/2, which 
implies that the upper bound on the temperature is T < T sta ti c - This results from the fact 
that in this case there is no geometric upper bound on k. Analogously to the S 5 case we find 

3k 

*±(r,k) = -__=, (6.9) 

a/1 - 8A±(a) 

with £l± given by Eq. (4.17) and r = r/L. 
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Thermodynamic properties 

In parallel with (5.1) the energy E, angular momentum J and entropy S for this case are 
computed to be 



E(r) = n {3) r 3 [ Rl 



e r 



Qr 



s ,3r2^(r) C;> . l _ n ,;; 



, J(r) = n (3) r 3 L 2 



, S(r)=n (3) r 3 s(r) 



(6.10) 



with e(r), s(r) given in (5.2). The resulting on-shell Gibbs free energy is F = — fi( 3 )(r 3 kP + 
Qr 4 /L), and by varying this for constant T, Q, and Q one can obtain the corresponding EOMs 
(6.7). Using the definitions of Eq. (5.5), we compute the following thermodynamic quantities 



E±(T,k) 



S±(T,k) 



where Cl = QL. 



n 



± 



J±(i + f|) 

k 5 f| 
4-^/5" T 5 cosh 4 a 



f 4 



J±(T,k) 



27 k 3 riOt 5 + 4sinh 2 a 



27 



, F±(T,k) 



" f 4 cosh 4 a 

27 k 5 f| 1 + 4 sinh 2 a 
16X/5 1 f 4 cosh 4 a 



(6.11) 
(6.12) 
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Figure 9: J versus f for the two solution branches of thermal giant gravitons on AdSs for 
f = 0,0.1,0.2 and 0.3. 

We have plotted the angular momentum J as function of the S s radius f in Fig. 9 as well 
as the free energy F as function of J in Fig. 10. It is clear from Fig. 9 that for T > the 
angular momentum is bounded from below as J > J m in > as in the case of the thermal 
giant graviton moving on S 5 . Correspondingly, the S 3 radius f is also bounded from below. 
Instead, there is no upper bound on J as in the extremal case. 



For J > J m i n there are two available solutions. From Fig. 10 we see that the stable 
solutions belong to the lower branch but the part of the lower branch between f j min and r m [ n 
is unstable. All solutions in the upper branch are unstable, as in the extremal case 12 . 

12 An off-shell analysis in the same spirit of Sec. 5 using the method described in App. D can be carried out 
giving further evidence for this. 
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0.2 0.4 0.6 0.8 1 

Figure 10: F versus J for the two solution branches of thermal giant gravitons on AdSs for 
T = 0, 0.1, 0.2 and 0.3. Note that the curves are very close to each other for different T. 

The analysis of the validity of the probe approximation works in the same way as for the 
S 5 thermal giant graviton. For small giant gravitons with r/L small we find the condition 
r s <C r, or equivalently, Nd^/N <C r 4 /L 4 . Instead when r/L is not small the condition is 
r s <C L which corresponds to A^D3 "C N. Note that the requirement T < Thp implies that 
T <C T stat i c . 

Low temperature expansion 

Focussing on the stable branch of the solution we find that for low temperatures 

n_(f,r) = ±- + 0(f 8 ) , 

J_(f,f) = r 2 + y^f 4 + 0(f 8 ), (6.13) 

F_(f,r) = f 2 + ^ 3 ^ 2 f 4 + 0(f 8 ) , 
27 r z 

Inverting the second equation and plugging into the third equation of (6.13) we find exactly 
the same result for F_(T, J) as given in (5.12). It is interesting that this leading thermal 
correction to the free energy in the low temperature limit is thus universal for both the lower 
branches of the thermal giant gravitons on S 5 and AdSs. 

7 Discussion and outlook 

In this paper we constructed and studied thermal giant gravitons. Thermal giant gravitons 
result from heating up the giant gravitons and the background space-time they move in. We 
focussed on giant gravitons that in the extremal case are 1 /2 BPS with < J < N (along with 
a companion non-BPS branch of solutions with N < J < 9N/8) obtained from D3-branes 
wrapped on an S 3 and moving in the AdSs x S 5 background. The thermal giant gravitons 
are described using black D3-brane probes as described by the blackfold approach [3, 1, 2]. 
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Using the AdS/CFT correspondence our thermal giant graviton solution is expected to 
correspond to a thermal state in the dual gauge theory. It would be highly interesting to find 
a description of this thermal state in the gauge theory and compare its properties to those of 
the thermal giant graviton. In particular, it would be important to compute the free energy 
correction on the gauge theory side that corresponds to our result Eq. (1.1) for the thermal 
giant graviton. 

An important question that arises from our work is what happens to the system of the 
thermal giant graviton in AdSs x 5 5 when the temperature is heated up beyond the Hawking- 
Page temperature Thp. For T < Thp the background is in the phase where thermal AdSs 
times 5 5 is dominant. Instead for T > Thp the dominant phase is an AdS black hole times S 5 . 
This means the giant graviton moves in the latter background for T > Thp. Thus, it would be 
interesting to repeat the analysis of this paper for the AdS black hole times S 5 background. 

One can also study what happens if one considers many giant gravitons moving in AdSs x 
S 5 . If sufficiently many giant gravitons move along the equator of 5 5 one can describe them 
as smeared along this circle. However, while for extremal giant gravitons the descriptions 
of the smeared and single giant gravitons are equivalent, for thermal giant gravitons the 
horizon topology would change as one increases the number of giant gravitons to the point 
where the horizons of each of them overlap. Thus, the non-BPS nature of thermal giant 
gravitons makes it particularly interesting to study the difference between the smeared and 
non-smeared phases. Moreover, the smeared phase is connected to the superstar [27] and 
1/2 BPS bubbling AdS solutions of LLM [28]. Indeed, we expect that the smeared thermal 
giant graviton solution should be a finite temperature version 13 of the 1/2 BPS bubbling AdS 
solution. 

Another direction to pursue would be to compute higher-order corrections 14 in the matched 
asymptotic expansion to our thermal giant gravitons. In this paper we have worked in the 
probe approximation Nn^ C iV to the leading order in the expansion parameter Nd^/N. 
Computing the higher-order corrections would reveal information on what happens for larger 
values of Nd^/N. This could possibly also be interesting to examine numerically. 

Finally, we have focussed in this paper solely on heating up 1/2 BPS giant gravitons. An 
interesting future direction to pursue would be to heat up giant gravitons with less supersym- 
metry [31]. Note that even in the extremal case not much work has been done to find the 
explicit brane configurations for 1/4 and 1/8 BPS giant gravitons. 
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A Detailed stability analysis of extremal giant gravitons 

In this appendix we carefully analyze the stability of the two solution branches of extremal 
giant gravitons by introducing a time dependent perturbation in the radial dynamics and solv- 
ing the linearized EOMs. More precisely, given a giant graviton configuration characterized 
by a size r andangular velocity O and with its dynamics governed by the DBI Langrangean 
(2.6), we consider the following perturbation around an on-shell configuration 

r = f (l + eC r e iuJT ) , Q = Q (l + eC n e luJT ) . (A.l) 

Here fo and f^o denote the on-shell values of the configuration we perturb around and the 
perturbation parameter e satisfies The factors C r and measure the amplitude of 

the oscillations. Our aim is to search for solutions with cu 2 > 0, signifying the stability of the 
configuration. 

In order to introduce the perturbations (A.l) in Eq. (2.6) we need to evaluate the quantities 
involved to lowest order in e. We begin by obtaining the induced metric 7 a 6 for this time- 
dependent embedding in the form 

jabda a da b = -k 2 dr 2 + r 2 dtt 2 {3) + 0(e 2 ) , (A.2) 

where, despite of their resemblance to the unperturbed case, k and r are now time-dependent 
quantities. Doing the same exercise for the WZ term of the action we obtain the time depen- 
dent Lagrangean 



Lr>Bl = — ?D3%3) 1 



2^.2 



L z r 



rQ 



(A.3) 



2p 2 k 

where we have defined p 2 = L 2 — r 2 . The two Euler equations for the Lagrangean above take 
the form 

7\2/ r 2 \ 

k 2 + r 2 ^ 2 - ArQk + rp^VLrti + rpk 2 r ((2 - L 2 Q 2 )r 2 - 3L 2 k 2 ) r 2 = , (A.4) 

p d k 2 \ 2 ' ) 

which describes the radial dynamics, while angular momentum conservation is encompassed 
by 

d ( 4 r 2 r 3 L 2 Q i p 2 r 3 Q 



In the EOMs (A.4), (A. 5), we introduce the perturbations (A.l) and we take the on-shell 
value Oq to be that of the lower or upper branch solution in (2.10). 
Starting with the lower branch, Eqs. (A.4) and (A. 5) reduce to 

{p 2 LC n f + 2r C r )e + O(e 2 ) = , Uc n + C I L ^j e + Q(e 2 ) = , (A.6) 

and are solved by ui 2 = 4/L 2 . Therefore, since u 2 > for all values of r, we conclude that the 
lower branch is always stable. Due to the 1/2-BPS property of the branch, this is expected. 
For upper branch we instead find the set of equations 

-2Cnf + 3 ^f-)e + O(e 2 ) = 0. (A.7) 
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(l8C n - C r r O (3LV + 16/rf2g)) e + C(e 2 ) = , (A.8) 

which are solved for 

Hence we conclude, as stated in Sec. 2.2, that the non-BPS upper branch is stable for r > 
\pSLj2 and unstable for r < 



B Thermodynamic blackfold action and Smarr relation 

In this appendix we show that the (mechanical) action (3.17) is equivalent to the thermody- 
namic action (3.18). To this end we first rewrite (3.17) as 

1 = At f dV {p) [C (hl) + £ {em) ] , (B.l) 

J Bp 

where from now on the subscripts "bf and "em" refer to the blackfold and external field 
respectively. In (B.l) we have factored out the integration over the (Killing) time t. This 
produces a redshift factor which must be included in the Lagrangian densities, e.g. £(bi) = 
■j^RqP where 7j_ is defined in (3.15). From the conserved quantities derived in (3.16) we 
also introduce the Hamiltonian and angular momentum densities 

H = n (M) +H {cm) = 7J 1 (t^ {) + Vj^j) n^u , J = J(bi)+J(em) = 7I 1 (rgf) + V ( 7 m) ) n^xu , 

(B.2) 

where is the blackfold stress tensor which encapsulates the gravitational and electro- 
magnetic self-energy /momentum and V^s (see (3.14)) is associated with the coupling of the 
charge current to the external electromagnetic field. Notice that the electromagnetic contri- 
butions only depend on the embedding degrees of freedom of the blackfold and not on the 
effective blackfold fluid degrees of freedom. 

Now, for the blackfold degrees of freedom we have the relation 

^(bf) + l^u^n^Ts = Q J (M) - £ (bf) , (B.3) 

which follows from eq.(2.19) of [6] by multiplying with 7j^ 1 . This is the blackfold generaliza- 
tion of the usual relation H = 9 J — L in Hamiltonian mechanics between Hamiltonian and 
Lagrangian, but now with an extra term contributing to the energy due to the fact that the 
blackfold has internal thermal degrees of freedom living on it. However, since the external 
electromagnetic field does not couple to the thermal degrees of freedom living on the blackfold, 
one has for the electromagnetic part that 

%(em) = &J(em) ~ £(em) ■ ( B - 4 ) 

We now use (B.3), (B.4) in (B.l) along with the expression (3.19) for the total entropy S of 
the blackfold. If we also rotate to Euclidean time so that At — > At = (5 = 1/T, we then find 
that the Euclidean action is given by 

I E = E-nj-TS. (B.5) 
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Here we recall that, as in [6] this is the Euclidean action at fixed charge Q p . As explained in 
that reference, it is also possible to go to an ensemble where the charge can vary by introducing 
a potential <£ p dual to the charge and performing a Legendre transformation. 



Smarr relation 

Finally, we derive the Smarr formula for blackfolds in external fields. We use the perfect fluid 
stress tensor T^ u = (e + P)u^u v + Ph^ and the local thermodynamic relations for charged 
p-branes inL> = n + p + 3 dimensions 

e + P = Ts, e = -(n + l)P -n^ p Q p . (B.6) 

First, we note that the Smarr relation found previously for blackfolds based on charged p- 
branes (with zero external field) is easily generalized to the case where £^ is not orthogonal 
to the worldvolume B p . One finds 

(D - 3)E [h{) — (D — 2) (OJ (bf) + TS) - n$> H Q p = , (B.7) 



where 



$H= f dV^^R^p , (B.8) 

J Br, 



7?b°f) = " / dV ( P ) hl 1 RotrT + J^T^nJ , trT = lab T ah . (B.9) 



B p 



We then add to both sides of (B.7) the term (D — 3)_EV em -) — (D — 2)QJ( em ), yielding the 
generalized Smarr relation 

(23 _ 3)2? — (D — 2) (OJ + TS) - n<5> H Q p = %* , (B.10) 

where 

T tot = - [ dy (p) ( 7 I%trr + 7 J 1 (r^ ) +V^ m p^ + (^-2)>C (em) ) . (B.ll) 

J Bp 

Note that, as expected, the total tension gets modified by the presence of the external field. 

C Derivation of the thermal giant graviton equation of motion 



In this appendix we provide some details on the derivation of the thermal giant graviton 
blackfold EOM (4.4). First we note that for a perfect fluid stress tensor, the extrinsic blackfold 
EOM (3.7) can be rewritten as [2, 6] 

PK" + sTii^ = T» , (C.l) 

where P is the pressure, s the entropy density, 7~ the local temperature, the extrinsic 
curvature vector, the fluid acceleration and the external force. 
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We proceed by computing the various terms. Using (3.8) along with the the background 
(2.2) and embedding (2.3) we compute 

c _ ^ Vi K c__m4 K c__wM K c__ r j l l (C2) 

^tt - j2 ' A ll ~ rL 2 ' ^22 - rL 2 ' ^33 ~ £ 2 ' 

where 

Hi = sj L 2 — r 2 , fj,2 = r sin 6 , \xj> = r cos , (C3) 

and we recall that £ is the direction parameterizing the fibration of the S 5 into S 3 s on which 
the giant graviton of size r = L sin £ is defined. This gives for the extrinsic curvature vector 

K C- Mi ^ 2 + 3k 2 

with the other components of equal to zero. Notice that is manifestly negative. To 
compute the extrinsic blackfold force term, we first compute the pull back of the RR field 
strength (2.1) on the blackfold according to ^123^ = 4£ld£l<$\/L. Using this together with 
the blackfold four-current (4.3) in the right side of (3.7), we find that the only non- vanishing 
component of the force term is 

* - ^ ■ (C5) 

We see that F^/(QQ) is manifestly positive. Finally we compute the fluid acceleration ii^ = 
u v V . To this end notice that although the (local boost) vector field u a is only defined on 
the world volume of the giant graviton, it can be pushed forward to the vector field on the 
entire RxS 5 . This means that the acceleration can be computed using it M = log k so that 

U = L2k2~ ' (C - 6) 

with the rest of the components of u M equal to zero. Note that the acceleration vS is manifestly 
positive. 

The extrinsic equation (C.l) is thus only non-trivial for [i in the ^-direction and from the 
results above we see that since the blackfold pressure P (see (4.1)) is negative, the left hand 
side is manifestly positive. The structure is therefore clear: In order for the D3 brane not to 
collapse under gravity, the electromagnetic repulsion term (the right hand side) must exactly 
balance the gravitational pull (the left hand side). We therefore conclude that QQ > 0, which 
means that the solution must always be rotating and charged, as expected. Moreover, we see 
that the solution space is symmetric under charge conjugation and time inversion. 

The thermal giant graviton blackfold EOM (4.4) of the text is then obtained from (C.l), 
by substituting (C.4), (C.5), (C.6) and using the pressure P in (4.1) along with T and s in 
(4.2). 



D Analysis of the two meeting points 

In this appendix we perform a careful analysis of the two meetings points in configuration 
space mentioned in Sec. 5 corresponding to the maximal and minimal charge parameter giant 
gravitons. These bifurcation points (see e.g. Ref. [26]) deserve special attention as they provide 
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key information about the overall stability properties of the thermal giant graviton. Being 
bifurcation points where two sets of equilibria configurations meet, the stability properties of 
the physical system can in generally change and therefore special attention is needed. 

According to the stability analysis that was presented in Sec. 5 for the low temperature 
regime, all configurations near the maximal size are stable in both branches. However, one 
can imagine increasing slightly the temperature and moving away from such regime while still 
being below the Hawking-Page temperature. The analysis of the stability properties of the 
maximal giant graviton carried out below allow us to conclude that, even after the temperature 
is further increased, the thermal gravitons are stable. Moreover, part of the lower branch, 
which contains J m i n , is not covered by the low temperature expansion. Since we do not reach 
J m i n analytically in that expansion, we instead analyze here the end point of the lower branch 
described by the minimal charge parameter giant graviton and show explicitly that a change 
of stability has occurred. We begin by describing these two different limiting cases and then 
proceed to study their stability in Sec. D.3. 

D.l The maximal giant graviton limit 

In this section we probe the regime in solution space dominated by the dynamics of the 
maximal giant graviton, which is achieved by performing an expansion of Eqs. (4.16)-(4.17) 
around k = 1. The extremal maximal giant graviton is the object for which there is good 
evidence that its dual operator is a Schur polynomial, hence this configuration is expected to 
be an ideal candidate for a preliminary study on the dual state of giant gravitons at finite 
temperature. In view of this, we begin by describing the properties of the thermal maximal 
giant graviton. At the exact point k = 1 (or r = L) the angular velocity of the configuration 
is given by 

where (fro is given in terms of 5q is obtained by setting k = 1 in Eq. (4.14), (4.15). The 
conserved charges can be obtained from Eqs. (5.1). The total energy and angular momentum 
then take the form 

E ± (0 O ) = ^^(4 + ^), J±(0o) = l, (D.2) 



where we see clearly that, as in the extremal case, the angular momentum is independent of 
the angular velocity and temperature of the configuration. Moreover the total entropy and 
Helmholtz free energy read 

s u \ - 27 $ p±u ^ - 27 ( 4 ~ 3 ^Q m v 

S±(0o) " igW^ ' Fh(0o) " i6v^ • (D ' 3) 



The large giant graviton expansion 

We now perform the expansion around k = 1. From Eq. (4.16) we find a relation between the 
values of f and the values of k of the form 

k- 1 , 
r± = 1 + . D.4 
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Parametrizing the expansion in terms of r we invert the above relation in order to find: 

(k±-l) = n|(0o)(f-l) . (D.5) 

Since the expansion is only valid for values of k ~ 1 we must require L 2 Q±((fto)(r - 1) < 1. 
Defining df = f — 1 and using Eq. (4.17) to perform the same expansion yields 



where we have defined the function g((f>o) through the expression: 

g(M = -8(1 + A±(0 O ))(-1 + A±(0o)) + 363 k A±(0 o ) . 
Similarly, the charge parameter <ft given in Eq. (4.15) is expanded to 

4> = <Po (l - &±(<Po)4>of(<Po)df^J , 
where we have defined the function f(4>o) as 

/(0o) = 2 (^2^/3 cos(^o) + VST cos(^o) + 8cos(^)sin 3 (^ sec(<5 )csc(c5o) 

The physical properties can be easily obtained from expressions (5.1) and read 

-3" 



E± = E± ( 



3 i + oiw -niw/w 



df 



4 + ^o 

J± = J±(0o) + (4 - 2ft±(0 o )0o(4 + O )J dr , 

S± = S±(0 O ) (l + [3 (l + fi|(0 o )) + 2^(0 O )(1 - f(M) 
and finally the Helmholtz free energy 

-12 + 90o + fi|(^o)(-12 + 170o - 54(2 - 30 o )/(0o)) 



dr 



F± = F±(0o)(l + 



4-30o 



(D.6) 

(D.7) 
(D.8) 

(D.9) 

(D.10) 

(D.ll) 
(D.12) 

df j .(D.13) 



D.2 The minimal charge parameter limit 

The minimal charge parameter limit is the point in configuration space where the two branches 
connect smoothly and can be seen as the point particle analog at finite temperature. Further- 
more, it is a limit which describes thermal giant graviton configurations at any temperature 
T, including very low temperatures, but it is not captured by the low temperature limit of 
Sec. 5 because f ~ T. At the exact meeting point the thermal giant graviton is characterized 
by a specific value of the charge parameter a which implies a characteristic size r given by 
(4.18) and corresponding angular velocity 



n 



1 + -T 2 . 

5 



The physical properties are then easily obtained from Eqs. (5.1) and read: 



E 



18 f 2 
5 & 



2P - r 4 



f 4 ^, F H 



9 



18 5 - 2T 2 , 



25 



h 2 



(D.14) 



(D.15) 
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The minimum charge parameter expansion 

We now obtain an effective description of the physics near this limit by expanding the physical 
properties of these configurations near k = T. We begin by expanding the charge parameter 
parametrized by <f> using Eq. (4.15), which to leading order yields 



9 -C I 75-1) , ( D -!6) 



where C is a numerical constant (~ v3 ) which can be calculated exactly. Using this form of 
the charge parameter we expand the size r using Eq. (4.16), which leads to 



r± = f\l + -^-= ^J^-l). (D.17) 



Thus we invert the equation above to find the relation 



II = ±lr *-, (D.18) 



where we have defined the expansion parameter 

^rbv 1 -?)- (D - 19) 

Therefore for this approximation to be valid we need to require df < 1. Using the above 
definition we can rewrite the correction to the angular velocity to leading order in the form 

- / 9T 2 A 
0± = n 1 ± — df . (D.20) 

V 5^ 2 / 

It is then straightforward to write down all the physical properties using Eqs. (5.1) 

e = e (i ± 3 ^,-i> d A , j = j(i± v.**** , (D . 21) 

V & J V n 2 (w -f 2 ) / 

S _ fi (l ± 6 <l^i (i A , F„ - F„ f 1 ± 3(5 + 2^-77") A 
V / V 2 ( 5 _ 2 f2) / 

D.3 Stability properties in the various limits 

In this section we examine the stability properties of the various limits of the thermal giant 
graviton taken in the previous sections. To this aim we consider the localized giant graviton 
to be in thermodynamical equilibrium with the surroundings at temperature T. Moreover, 
since the total angular momentum J is conserved, the relevant variables for describing the 
thermodynamic ensemble are the size of the giant graviton r, the temperature T, the angular 
momentum J and the (conserved) total charge Q = Id3-/Vd3- The stable solutions to the 
blackfold EOMs are then characterized by the paths in configuration space for which the 
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Helmholtz free energy Fn = E — TS is minimized for T, J and Q held fixed. In other words, 
the stable solutions are determined by the requirements 



( F h)(i)-^_ = and (F H ) (2) ^^_>0. (D.23) 



t,j,q ' ' K ' dr 



TJX 



The first of these equations is equivalent to the EOM (4.4) and was examined in Sec. 4. The 
formulae (5.1) for the conserved quantities allows us to obtain the free energy of a (in general 
off-shell) given thermodynamical configuration. However, notice that the free energy will 
be parameterized in terms of the angular velocity f2 and charge parameter <j). Determining 
the derivatives (D.23) is straight forward. For a given r, T and Q, let = Q(r;T,Q) and 
(p = (f)(r; T, Q) denote the corresponding on-shell values. Now consider a small variation of the 
configuration r — > r + dr, Q — > 0(r; T, Q) + fi( 1 )dr + f2( 2 )dr 2 , cj) — > 4>{r; T, Q) +<j)^dr + 0( 2 )dr 2 
so that J, T and Q are kept constant up to 0(dr 2 ). This gives us four equations (using 
equations (5.1) and (4.6) for each order) which allows us to determine the four parameters 
JVn, fi(2)> <^(i)) <f>(2)- Since we are perturbing around an on-shell configuration (Fh)(i) = 0, so 
that the overall change in Fr is 

F H ^F U + (F H ) (2) dr 2 . (D.24) 

Inspecting the sign of (i ? H)( 2 ) allows us to determine the stability of the given solution to the 
blackfold EOM. This method was used in Sec. 5.3 to probe the stability properties in the low 
temperature regime. We now make further use of this method by applying it to the two cases 
presented in the previous sections and state the results for the second order change in the free 
energy in the various limits. 

Maximal giant graviton limit 

Close to maximality we make use of the expansion given in Sec. D.l. We find for the second 
order change in the free energy 



v± 4 4-4 ( Ao-0 ± (^o)^^W(4 + ^o) + ^(^o)(l + f) s 
(F H ) (2) = ^ n ± '(a i 1 \ + H (<fo> r ) dr ' ( D - 25 



where H{4>o,r) is some intricate function of 4>q and f. An explicit computational check can 
be performed in order to conclude that (Fh)^ 2 ) > f° r an values of r except at T = 1 (see 
below). Hence, near maximality, all giant graviton configurations are stable and so remain 
for all temperatures. 

Minimal charge parameter limit 

At the bifurcation point in configuration space where the two branches meet smoothly, the 
expansion given in Sec. D.2 is the required tool to study the stability properties of these 
configurations. Using the same method as for the previous case we find the second order 
change in the free energy: 

(F„)f 2) = (F„) (2) (l ± , (F H)(2) = -^F H . (D.26) 
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Careful inspection of the sign of (Fh)^ leads us to conclude that all giant graviton configura- 
tions near the bifurcation point are unstable, i.e., (Fh)( 2 ) < 0- This provides further evidence 
for the change of stability properties occurring at J m i n and J max as explained in Sec. 5. 
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